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Abstract —Heterogeneous networks (HetNets) have strong in¬ 
terference due to spectrum reuse. This affects the signal-to- 
interference ratio (SIR) of each user, and hence is one of the 
limiting factors of network performance. However, in previous 
works, interference management approaches in HetNets are 
mainly based on interference level, and thus cannot effectively 
utilize the limited resource to improve network performance. 
In this paper, we propose a user-centric interference nulling 
(IN) scheme in downlink two-tier HetNets to improve network 
performance by improving each user’s SIR. This scheme has 
three design parameters: the maximum degree of freedom for 
IN (i.e., maximum IN DoF), and the IN thresholds for the macro 
and pico users, respectively. Using tools from stochastic geometry, 
we first obtain a tractable expression of the coverage (equiva¬ 
lently outage) probability. Then, we characterize the asymptotic 
behavior of the outage probability in the high reliability regime. 
The asymptotic results show that the maximum IN DoF can 
affect the order gain of the asymptotic outage probability, while 
the IN thresholds only affect the coefficient of the asymptotic 
outage probability. Moreover, we show that the IN scheme can 
linearly improve the outage performance, and characterize the 
optimal maximum IN DoF which minimizes the asymptotic 
outage probability. 

1. Introduction 

The modern wireless networks have seen a significant 
growth of high data rate applications. The conventional cel¬ 
lular solution, which comprises of high power base stations 
(BSs), cannot scale with the increasing data rate demand. 
One solution is the deployment of low power small cell BSs 
overlaid with conventional large power macro-BSs, so called 
heterogeneous networks (HetNets). HetNets are capable of ag¬ 
gressively reusing existing spectrum assets to support high data 
rate applications. However, spectrum reuse in HetNets causes 
strong interference. This affects the signal-to-interference ratio 
(SIR) of each user, and hence is one of the limiting factors of 
network performance. Interference management techniques are 
thus desirable in HetNets. One such technique is interference 
cooperation. For example, in [1], a cooperation strategy among 
a fixed number of strongest BSs of each user is investigated for 
HetNets. Reference [2] proposes another interference cooper¬ 
ation strategy among BSs whose long-term average received 
powers (RP) are above a threshold. However, both cooperation 
strategies do not take into account the strength of each user’s 
desired signal, and thus cannot effectively improve each user’s 
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SIR. Moreover, both [1] and [2] only consider single-antenna 
BSs. 

Deploying multiple antennas at each BS in HetNets can 
further improve network performance. With multiple antennas, 
besides boosting the desired signal to each user, more effective 
interference management techniques can be implemented. In¬ 
terference coordination strategies using multiple antennas have 
been extensively investigated in single-tier cellular networks 
(see [3,4] and the references therein). However, there are 
very limited works on the analysis of interference manage¬ 
ment techniques in multi-antenna HetNets. In [5], the authors 
consider a HetNet with a single multi-antenna macro-BS 
and multiple small-BSs, where the multiple antennas at the 
macro-BS are used for serving its scheduled users as well 
as mitigating the interference to the receivers in the small 
cells using different interference coordination schemes. These 
schemes are analyzed and shown to improve the performance 
of the HetNet. However, since only one macro-BS is con¬ 
sidered, the analytical results obtained in [5] cannot reflect 
the macro-tier interference, and thus cannot offer accurate 
insights for practical HetNets. In [6], a fixed number of BSs 
form a cluster, and adopt an interference coordination scheme 
in downlink multi-antenna HetNets, where the BSs mitigate 
interference to users in the cluster. Bounds of the coverage 
probability are derived based on an unrealistic assumption 
that the BSs in each cluster are the strongest BSs of all the 
users in this cluster. This scheme cannot effectively improve 
each user’s SIR, as it does not consider each user’s desired 
signal strength. Recently, a novel user-centric IN scheme, 
which takes account of each user’s desired signal strength 
and interference level, is proposed and analyzed for single¬ 
tier cellular networks with multi-antenna BSs in [4]. However, 
in [4], the maximum degree of freedom for IN (IN DoF) at 
each BS is not adjustable. This scheme thus cannot properly 
exploit the limited resource in single-tier cellular networks. 
Moreover, directly applying the approach in [4] to HetNets 
cannot fully exploit different properties of the macro and pico 
users in HetNets. 

In this paper, we consider the downlink two-tier multi¬ 
antenna HetNets and propose a user-centric IN scheme to 
improve network performance by improving each user’s SIR. 
In this scheme, each macro (pico) user first sends an IN request 



to a macro-BSj^if the power ratio of its desired signal and the 
interference from the macro-BS, referred to as the signal-to- 
individual-interference ratio (SIIR), is below an IN threshold 
for macro (pico) users. Then, each macro-BS utilizes zero¬ 
forcing beamforming (ZFBF) precoder to avoid interference 
to at most U users which send IN requests to it as well as 
boost the desired signal to its scheduled user. This scheme has 
three design parameters: the maximum IN DoF U, and the IN 
thresholds for the macro and pico users, respectively. In gen¬ 
eral, the investigation of interference management techniques 
in multi-antenna HetNets is very challenging, mainly due to 
i) the statistical dependence among macro-BSs and pico-BSs 
[5], ii) the complex distribution of desired signal using multi¬ 
antenna communication schemes, and iii) the complicated 
interference distribution caused by interference management 
techniques (e.g., beamforming). In this paper, by adopting 
appropriate approximations and utilizing tools from stochastic 
geometry, we first obtain a tractable expression of the coverage 
(equivalently outage) probability. Then, we characterize the 
asymptotic behavior of the outage probability in the high reli¬ 
ability regime. The asymptotic results show that the maximum 
IN DoF can affect the order gain of the asymptotic outage 
probability, while the IN thresholds only affect the coefficient 
of the asymptotic outage probability. Moreover, we show that 
the IN scheme can linearly improve the outage performance, 
and characterize the optimal maximum IN DoF which mini¬ 
mizes the asymptotic outage probability. The analytical results 
obtained in this paper provide valuable design insights for 
practical HetNets. 

II. System Model 

A. Downlink Two-Tier Heterogeneous Networks 

We consider a two-tier HetNet where a macro-cell tier is 
overlaid with a pico-cell tier, as shown in Fig.[2 The locations 
of the macro-BSs and the pico-BSs are spatially distributed 
as two independent homogeneous Poisson point processes 
(PPPs) and <I >2 with densities Ai and A 2 , respectively. The 
locations of the users are also distributed as an independent 
homogeneous PPP with density A^^. Without loss of gener¬ 
ality, denote the macro-cell tier as the 1st tier and the pico-cell 
tier as the 2nd tier. We focus on the downlink scenario. Each 
macro-BS has Ni antennas with total transmission power Pi, 
each pico-BS has N 2 antennas with total transmission power 
P 2 , and each user has a single antenna. We assume > N 2 . 
We consider both large-scale fading and small-scale fading. 
Specifically, due to large-scale fading, transmitted signals from 
the jth tier with distance r are attenuated by a factor 
where aj > 2 is the path loss exponent of the jth tier and 
j = 1,2. For small-scale fading, we assume Rayleigh fading 
channels. 

^Note that, compared to pico-BSs, macro-BSs usually cause stronger 
interference due to larger transmit power, and each macro-BS normally has 
a larger number of transmit antennas (i.e., a better capability of performing 
spatial cancellation). Thus, it is more advisable to perform IN at macro-BSs. 


B. User Association 

We assume open access [1]. User i (denoted as Ui) is 
associated with the BS which provides the maximum long¬ 
term average RP among all the macro-BSs and pico-BSs. 
This associated BS is called the serving BS of user i. Note 
that within each tier, the nearest BS to user i provides the 
strongest long-term average RP in this tier. User i is thus 
associated with (the nearest BS in) the j*th tier, i^ j* = 
arg maxj^i^ where Zij is the distance between 

user i and its nearest BS in the jth tier. We refer to the 
users associated with the macro-cell tier as the macro-users, 
denoted as Ui = > P 2 Z^ 2 ^^, and the users 

associated with the pico-cell tier as the pico-users, denoted 
as U 2 = {ui\P 2 Z^ 2 ^ > PiZ^^^^Y All the users can be 
partitioned into two disjoint user sets: Ui and U 2 . After user 
association, each BS schedules its associated users according 
to TDMA, i.e., scheduling one user in each time slot. Hence, 
there is no intra-cell interference. 

C. Performance Metric 

In this paper, we study the performance of the typical 
user denoted as uq, which is a scheduled user located at the 
origin [2]. Since HetNets are interference-limited, we ignore 
the thermal noise in the analysis of this paper. Note that the 
analytical results with thermal noise can be calculated in a 
similar way. We investigate the coverage probability of uq, 
which is defined as the probability that the SIR of uq is larger 
than a threshold [1], and can be mathematically written as 

5(/3) = Pr (SIRo > /3) (1) 

where [3 is the SIR threshold. 

HI. Interference Nulling Scheme 

In this section, we first elaborate on a user-centric IN 
scheme to avoid interference from some macro-BSs which 
generate dominant interference. Then, we obtain some dis¬ 
tributions related to this scheme. 

A. IN Scheme Description 

First, we refer to a macro-BS which causes the SIIR at a 
scheduled user i in the jth tier {ui G Uj) below threshold Tj > 
1 as a potential IN macro-BS of Ui, where j = 1,2. We refer to 
Tj as the IN threshold for the jth tier. Mathematically, macro- 
BS ^ is a potential IN macro-BS of scheduled user Ui G Uj if 

Pi\Di ^ distance from macro- 

BS ^ to Note that Ti and T 2 are two design parameters of 
the IN scheme. In each time slot, each scheduled user sends IN 
requests to all of its potential IN macro-BSs. We refer to the 
scheduled users which send IN requests to macro-BS I as the 
potential IN users of macro-BS I (in this time slot). Consider a 
particular time slot. Let Ki denote the number of the potential 
IN users of macro-BS L Consider two cases in the following, i) 
If > 0, macro-BS I makes use of at most U (0 < U < A^i) 

^In the user association procedure, the first antenna is normally used to 
transmit signal (using the total transmission power of each BS) for RP 
determination according to LTE standards. 
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between pico-BS I and its scheduled user. Note that the non- 
IN scheme can be included in the IN scheme as a special case 
by letting Ti = T 2 = 1 and U = o|^ 

We now obtain the SIR of the typical user. Under the above 
IN scheme, uq experiences three types of interference: 1) 
residual aggregated interference from its potential IN 

macro-BSs ^ic which do not select for IN, 2) aggregated 
interference Ij^io from interfering macro-BSs ^10 which are 
not its potential IN macro-BSs, and 3) aggregated interference 
Ij ^2 from all interfering pico-BSs. Specifically, the SIR of the 
typical user uq ^Uj (j = 1, 2) is given 


SIRj,o = 






PlIj,lC + Pllj,10 + P2lj,2 


Fig. 1. System Model {U = 1) 


DoF to perform IN to some of its potential IN users. Note that 
U is another design parameter of this IN scheme. In particular, 
if 0 < < U, macro-BS I can perform IN to all of its 

potential IN users using Ki DoF; if Ki > U, macro-BS £ 
randomly selects U out of its Ki potential IN users according 
to the uniform distribution, and perform IN to the selected U 
users using U DoF. Hence, if Ki > 0, macro-BS i performs 
IN to = min (U^Ki) potential IN users (referred to as 
the IN users of macro-BS i) using DoF (referred to as 
the IN DoF), ii) If K£ = 0, macro-BS i does not perform IN. 
In this case, we let = 0. In both cases, Ni — DoF 
at macro-BS £ is used for boosting the desired signal to its 
scheduled user. 

Now, we introduce the precoding vectors at macro- 
BSs in the IN scheme. Consider two cases in the fol¬ 
lowing. i) If Ki > 0, macro-BS £ utilizes the low- 

complexity ZFBF precoder to serve its scheduled user and 
simultaneously perform IN to its 14 in,o IN users. Specif- 
ically, denote = [hi,^ wher^ 

hi^i ^ CA/’ati,! (Otvixi, Iati) denotes the channel vector 
between macro-BS £ and its scheduled user, and ^ 

CJV ATi, 1 ( 07 Vi X 15 Iati ) denotes the channel vector between 
macro-BS £ and its ith IN user (i = 1,..., The 

ZFBF precoding matrix at macro-BS £ is designed to be 

^ and the ZFBF vector at macro- 

BS ^ is designed to be fi^i = where is the 

first column of Wi^^. ii) If K£ = 0, macro-BS £ uses 
the maximal ratio transmission (MRT) precoder to serve its 
scheduled user, which is a special case of the ZFBF precoder 
introduced for > 0 and can be readily obtained from it 
by letting = 0, i.e., = h\ Next, we introduce 

the precoding vectors at pico-BSs. Each pico-BS utilizes the 
MRT precoder to serve its scheduled user. Specifically, the 
beamforming vector at pico-BS £ is f 2 / = where 

f^ 2 ,i ^ CA/'Ar 2 ,i (Oatsxi? Iats) denotes the channel vector 
^The notation X i> Y means that X is distributed as Y. 


where hy^oo is the channel vector between uq G Uj 
and its serving BS Bj,o, Yj is the distance between uq 
and fj^o is the beamforming vector at Bj^q, with 


^i,ooU,o 


Gamma (My, 1), Mi = Ni — 14in,o and 


M 2 = N 2 . Here, /y,ic = 




P 


j^o 


E 


eei‘2\Bj,Q 


Efe4>io 

\D2,eo\~"^ 




1/0 I 
2 




and Ij ^2 = 


4/0^2/ 


, where hy^^o is the channel 


vector between BS ^ in the jth tier and uq, fj^£ is the 


beamforming vector at BS ^ in the jth tier, and 
Gamma(l, 1). 




B. Preliminary Results 


In this part, we evaluate some distributions related to the 
IN scheme, which will be used to calculate the coverage 
probability in Q- We first calculate the probability mass 
function (p.m.f.) of the number of the potential IN users of 
iio’s serving macro-BS (when uq denoted as Ffo- The 

p.m.f. of Ffo depends on the point processes formed by the 
scheduled macro and pico users, which are related to but 
not PPPs 1?]. For analytical tractability, we approximate the 
scheduled macro and pico users as two independent PPPs with 
densities Ai and A 2 , respectively]^ Then, we have the p.m.f. 
of Kq as follows: 

Lemma I: The p.m.f. of Kq is given by 

~r ^ 

Pr(Ko=fc)«-^exp(-Z) (3) 


where L = Li -\- L 2 with 


Lj =27rXj / r / , ^ fYj{y)dydr 

Jo 


i = l,2. (4) 


^All the analytical results in this paper hold for Ti = T 2 = 1 and U = 0. 
^The index of the typical user and its serving BS is 0. 

^Note that approximating the scheduled users as a homogeneous PPP has 
been consi dered in existing papers (see e.g., [7]). Moreover, simulation results 
in Section [TV-A| will verify the accuracy of this approximation. 


























Here, fy-iy) (j = 1,2) are given as follows: 


fvAv) = exp (-TT ( Ai 2 /^ + A 2 y 


fv^iy) = ^^ 2 /exp f-TT ( Ai y + \ 2 y‘^ 


where Aj = Pr (r^o ^ ^j) 0 = 2) are given in ([^ and 

at the top of the next page. 

Proof: See Appendix [A[ ■ 

Next, we calculate the p.m.f. of rtiN,o = min {U, Kq) based 
on Lemma [ 7 ] which is given as follows: 

Lemma 2: The p.m.f. of 14 in,o is given by 


Pr ('?^iN,o = u) = 


Pr {Kq = u) , for 0 < u < U 

EZu^^i^o = k) , for u = U 
Let Pc{U,Ti,T 2 ) denote the probability that a randomly 
selected (according to the uniform distribution) potential IN 
macro-BS of uq selects uq for IN. Note that the event 
that Uq sends IN requests and the event that all the other 
scheduled users send IN requests are dependent. For analytical 
tractability, we approximate these two events as independent 
events. Then, we have Pc{U,Ti,T 2 ) ~ E 
which can be calculated as follows: 

Lemma 3: The IN probability is given by 


'{t&i.'} 


/U—1 Yk ^ fk \ 

P.(Kr.,r.)=«p(-i) ^- + (7^— . 

/c—0 k—U / 

Proof: See Appendix]^ ■ 

Note that different potential IN macro-BSs of uq selects 
Uq for IN dependently (as the numbers of the potential IN 
users of these macro-BSs depend on the locations of these 
macro-BSs and are thus dependent). For analytical tractability, 
we assume that different potential IN macro-BSs of uq select 
Uq for IN independently. Using independent thinning, ^ic 
can be approximated by a homogeneous PPP with density 
P-, {U, Ti, T 2 ) Ai, where pg {U, Tj, T 2 ) = 1 - Pc {U, Tj, T 2 ). 


IV. Coverage Probability Analysis 

In this section, we investigate the coverage probability in 
the general and small SIR threshold regimes, respectively. 


A. General SIR Threshold Regime 


As discussed in Section |II-B[ the typical user uq is in one 
of two disjoint user sets: Ui and U 2 . By utilizing tools from 
stochastic geometry and the preliminary results obtained in 
Section III-B| we have the following theorem: 

Theorem 1 (Coverage Probability): Under design parame¬ 
ters U, Ti and T 2 , we have 


1) coverage probability of a macro-user: 
5i(/3,f/,Ti,T2) = Pr(SIRo >/?|mo eWi), given 

in (|^ at the top of the next page. 


2) coverage probability of a pico-user: ^2 (/5, U, Ti, T 2 ) = 
Pr (SIRo > P\uq G ^^ 2 ), given in at the top of the 
next page, 

3) overall coverage probability 5 (/3, U, Ti, T 2 ) = 

AiSi{i3,U,TyT2) + A2S2{l3,U,TyT2), where Ai 
and A 2 are given in 0 and ^ at the top of the next 
page. 


Here, cP (s, rj,io) and (s.rj^k) (k e {10,2}) 


(n) 




j,k 


are given in O. and at the top of the next page, 
respectively, (s,and Ci.ys,rj^k) 

are given in I® and •El at the top of the next page 
(with J(10) = 1 and J(2) = 2), respectively. More¬ 
over, B {a,b,z) = f^u^~^{l — u)^~^du (0 < z < 1), 
Vn = {(na)^=i|na € = ”•}’ Mn = 

{(ma)a=l|TOa ^ , Y3a=i a ■ nia = u}. 

Proof: See Appendix ■ 

Theorem allows us to easily evaluate the coverage prob¬ 
ability in a numerical way. Fig. 2(a) plots the coverage 
probability versus the IN thresholds Ti and T 2 . We see from 
Fig. |2(a)| that the “Analytical” curves, which are plotted using 
5 (/3, U, Ti, T 2 ) in Theorem^ are reasonably close to the 
“Monte Carlo” curves (the error is no larger than 5.3%), 
although The orem\I\ is obtained based on some approximations 
(cf. Section |III-B| r Compared to the non-IN scheme, we 
observe that the IN scheme can achieve a good coverage 
probability gain (up to 18.6%). Moreover, we see that the 
coverage probability depends on both Ti and T 2 . Specifically, 
when Ti (T 2 ) is small, increasing Ti (T 2 ) increases the chance 
of macro-BSs performing IN to near users; while when Ti 
(T 2 ) is large, increasing Ti (T 2 ) reduces the chance of macro- 
BSs performing IN to near users, as the resource is wasted 
in performing IN to users far away. Further, we see that the 
coverage probability can be improved by jointly adjusting Ti 
and T 2 . 


B. Small SIR Threshold Regime 


To obtain more insights on the impact of the design param¬ 
eters, in this part, we investigate the asymptotic behavior of 
the IN scheme in the high reliability regime where /3 ^ o[^ 

I) Asymptotic Outage Probability Analysis: In this part, we 
analyze the asymptotic outage probability when /3 ^ 0. First, 
we define the order gain of the outage probability [8]: 


logPr(SIRo<^) 

d = hm -^--- 

/3^0 log P 


( 15 ) 


Recently, a tractable approach has been proposed in [9] to 
characterize the order gain for a class of communication 
schemes in wireless networks which satisfy certain conditions. 
However, this approach does not provide tractable analytical 
expressions for the coefficient of the asymptotic outage prob- 


Cl {s,r) denotes the Laplace transform of the aggregated interference 1. 
^Note that small ^ is a practical regime. Lor example, in current standards 
(e.g., 3GPP LTE), /3 is very small (e.g., -7 dB) [8]. In addition, in wide-hand 
systems (e.g., CDMA and UWB), [3 is usually small [8]. 
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ability, i.e., 

Pr(SIRo <E) 

/3™ [3^ 


(16) 


for most of the schemes using multiple antennas in the 
class. By utilizing series expansion of some special functions 
and dominated convergence theorem, we characterize both 
the order gain and the coefficient of the asymptotic outage 
probability of the IN scheme in multi-antenna HetNets, which 
are presented as follows. 

Theorem 2 (Asymptotic Outage Probability): Under design 
parameters U, Ti and T 2 , when /3 ^ 0 , we hav^ 

1 ) outage probability of a macro-user: 1 — 

5i(/3,(/,Ti,T2) -61 ((/,Ti,T2)/3^i-^ 

2 ) outage probability of a pico-user: 1 — S2{dj (7, Ti, T 2 ) ~ 

b2iU,Ti,T2)d^\ 


3) overall outage probability: 1 — S {d,U,Ti,T 2 ) ^ 

b {U, Ti, T 2 ) /3min{iVi-;7,iV2}^ where 

b{U,TuT2) 

b2{U,Ti,T2), U<N,-N2 

= { 6i(C/,Ti,r2) + 62 (C/,Ti,T 2 ), U = Ni-N2. 

6 i([/,Ti,T 2 ), U>Nj-N2 


Here, bj ([/, Ti, T 2 ) is given in (19 ( at the top of the next page 
with 


Uj = 


U, ifj = l,Ti>l,T 2 >l 
0 , otherwise 


and 


Pj = 


Pr (uiN,o = U), if j = 1, Ti > 1, T 2 > 1 
1 , otherwise 


(17) 


( 18 ) 


^f{x) ^ 9 {x) means lima^^o = 1- 


Moreover, 62 (U, Ti,T 2 ) decreases with U. 




























(a) = 9, /3 = 10 dB 


the result of the order gain in Theorem extends the existing 
order gain result in single-tier cellular networks where the DoF 
for IN in the ZFBF precoder is deterministic [3]. 


Fig. |2(b)| plots the outage probability versus SIR threshold p 
for different design parameters. From Fig. |2(b)| we clearly see 
that the outage probability curves with the same maximum IN 
DoF U have the same slope (indicating the same order gain). 
For the two outage probability curves with the same U but 
different Ti and T 2 , we observe a shift between these two 
curves (indicating different coefficients). Therefore, Fig. 2(b) 
verifies Theorem and shows that the result in Theorem f2 
can effectively refiect the outage probability for small (3. 



(b) 

Fig. 2. Coverage and outage Probability at Ni = 10, N 2 = 8, = 4.5, 

Oi 2 = 4.7, = 15 dB, Ai = 0.0005 nodes/m^, and A 2 = 0.001 nodes/m^. 

Proof: See Appendix [P] ■ 

Note that the IN scheme has three design parameters: the 
maximum IN DoF (i.e., U) and the IN thresholds (i.e., Ti and 
T 2 ). From Theorem we clearly see that the maximum IN 
DoF and the IN thresholds affect the asymptotic behavior of 
the outage probability in dramatically different ways. Specifi¬ 
cally, the maximum IN DoF U can affect the order gain, while 
the IN thresholds can only affect the coefficient. In addition, 
we see that U affects the order gain of the asymptotic outage 
probability through affecting the order gain of the asymptotic 
macro-user outage probability. Note that in this paper, IN is 
only performed at macro-BSs, and U is the upper bound of 
the actual DoF for IN in the ZFBF precoder (which is random 
due to the randomness of the network topology). Therefore, 


2 ) Optimization of Maximum IN DoF: From Theorem 
we know that U has a larger impact on the asymptotic outage 
probability than the IN thresholds. In this part, we characterize 
the optimal maximum IN DoF which maximizes the order 
gain, and the optimal maximum IN DoF which minimizes the 
asymptotic outage probability, respectively. 

From Theorem we see that the order gain of the asymp¬ 
totic outage probability is min{A^i — f/, A^ 2 }- Thus, we also 
denote it as d{U). First, we introduce the optimal design 
parameter which maximizes d{U), i.e., 

= arg max d(U) 

= arg max minjA^i — U,N 2 }. (20) 

Thus, the optimal order gain of the asymptotic outage proba¬ 
bility is d{U^). Next, we introduce the optimal design parame¬ 
ter f/*(/3, Ti, T 2 ) which minimizes the outage probability, i.e., 

[/*(/3,Ti,T2)^arg min (1 - 5 Ti, T 2 )) . 

( 21 ) 

The properties of and f/*(/3, Ti, T 2 ) are given by: 
Lemma 4 (Optimality Properties): 


1) When /3 ^ 0, we have d{Uf) = N 2 for all ^ 

{0,1,...,7Vi-7V2}; 

















2 ) 3/3 > 0 , such that for all /3 < /3, we havj^ 


f/*(Ti,T2) 




N 1 -N 2 -I, ifb2{Ni-N2-l,Ti,T2) < 
hi {Ni - 7V2, Ti, T2) + 62 {Ni - 7V2, Ti, T2) . 
Ni — N 2 , otherwise 


( 22 ) 

Proof: See Appendix]^ ■ 

Result 1) in shows that is independent of U. 

Thus, the IN scheme does not provide order-wise performance 
improvement compared to the non-IN scheme. Note that 
the impact of on the asymptotic outage probability is 
through the coefficient. Result 2) of Lemma tells us that the 
asymptotic outage probability 1 — 5 (^5, T 2 ), Ti, T 2 ), 

optimized over U for given Ti and T 2 , can be further opti¬ 
mized by optimizing Ti and T 2 , i.e., 


min 1 

Ti>1,T2>1 


S{(i,U*{TuT2),T^,T2) 


(23) 


We shall investigate this optimization problem in future work. 


V. Conclusions 

In this paper, we propose a user-centric IN scheme in 
downlink two-tier multi-antenna HetNets. Using tools from 
stochastic geometry, we first obtaine a tractable expression of 
the coverage probability. Then, we characterize the asymptotic 
behavior of the outage probability in the high reliability 
regime. The asymptotic results show that the maximum IN 
DoF can affect the order gain of the asymptotic outage 
probability, while the IN thresholds only affect the coefficient 
of the asymptotic outage probability. Moreover, we show that 
the IN scheme can linearly improve the outage performance, 
and characterize the optimal maximum IN DoF which mini¬ 
mizes the asymptotic outage probability. The analytical results 
obtained in this paper provide valuable design insights for 
practical HetNets. 
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Appendix 


A. Proof of Lemma ^ 


According to Slivnyak’s theorem [10], we focus on a macro- 
BS located at origin. We refer to this macro-BS as macro- 
BS 0. Note that both scheduled macro-users and scheduled 
pico-users may send IN requests to macro-BS 0. We first 
characterize the probability that a scheduled macro-user sends 
an IN request to macro-BS 0. Denote Ru as the distance 
between macro-BS 0 and a randomly selected (according to 
the uniform distribution) scheduled macro-user which we refer 
to as scheduled macro-user i. According to Section |III-A 


scheduled macro-user i sends an IN request to macro-BS 0 
if 


PiY,- 


< Ti and < Ri 


(24) 


Hence, scheduled macro-user i sends an IN request to macro 
BS 0 with probability pii{Ti) which is given by 


Pii{Ti) = Pr ( Ru <Yi < Ru 


(25) 


Assume that the scheduled macro-users form a homogeneous 
PPP with density Ai. Conditioned on Ru = r, we have 
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^^As U*(^,Ti,T 2 ) is independent of /3, we also write it as U*(Ti,T 2 ). 


Pii,Ru{r,T^) = Pr {ti r<Y,<r^ 

= /-J- /yi(y)dy (26) 

where fviiy) is the probability density function (p.d.f.) of 
Yi given by < 0 . Then, the scheduled macro-user density 
at distance r away from macro-BS 0 is (r, Ti)Ai. 

This indicates that the scheduled macro-users which send IN 
requests to macro-BS 0 form an inhomogeneous PPP with 
density Ti)Ai at distance r away from macro-BS 0. 

Next, we characterize the probability that a scheduled pico- 
user sends an IN request to macro-BS 0. Denote R 2 i as 
the distance between macro-BS 0 and a randomly selected 
(according to the uniform distribution) scheduled pico-user 
which refer to as scheduled pico-user i. Similarly, we assume 
that the scheduled pico-users form a homogeneous PPP with 
density A 2 , and it is independent of the PPP formed by the 
scheduled macro-users. Then, we can show that the scheduled 
pico-users which send IN requests to macro-BS 0 form an 
inhomogeneous PPP with density (^5 ^ 2 )A 2 distance 





r away from macro-BS 0, where 


TABLE I 

Parameter values 


P 2 i,R 2 i{r,T 2 ) = Pr 


I 




PtT2 

L o 

2 j^OL2 


r «2 < I 2 < 


fYAv)<^y 


^2 r^2 


M ^2_V 

where fY 2 (y) is the p.d.f. of Y 2 given by 



According to the superposition property of PPPs [10], the 
scheduled macro-users and the scheduled pico-users which 
send IN requests to macro-BS 0, i.e., the potential IN users of 
macro-BS 0, still form a PPP with density (r, Ti)Ai + 

^ 2)^2 at distance r away from macro-BS 0. There¬ 
fore, the number of the potential IN users of macro-BS 0 is 
Poisson distributed with the following parameter (i.e., mean) 

POO 

L = 27r r (r, Ti)Ai + P 2 i,R 2 i (a ^ 2 )A 2 ) dr 

_ _ 

= Ii + I 2 . (28) 


B. Proof of Lemma ^ 


From Section |III-A| we know that whether a scheduled 
user sends an IN request to a macro-BS or not depends on 
its location relative to this macro-BS. Hence, the event that 
uq sends an IN request to one of its potential IN macro- 
BSs and the event that any other scheduled user sends an IN 
request to the same macro-BS are dependent. The dependence 
is especially high when uq and that scheduled user are close. 
For analytical tractability, we approximate these two events as 
independent events. Then, we have 


Pc {U, 


U-1 




= ^Pr(/fo = fc)+E 


u 


/c +1 


Fi{Ko = k) . (29) 


k=0 k=U 

Substituting <0 into ( [ 2 ^ , we have the final result. 


C. Proof of Theorem^ 

Let Rj^ic and Rj^io denote the minimum and maximum 
possible distances between uq G Uj and its nearest and furthest 
macro-interferers (among i^o’s potential IN macro-BSs which 
do not select uq for IN), respectively. Let Rj ^2 denote the 
minimum possible distance between uq G U 2 and its nearest 
pico-interferer. The relationships between Rj^ic^ 
and Yj, respectively, are shown in Table Based on 0 and 


j 

r.pic 

^.7,10 

^.7,2 

1 

Yi 

1 

T“i Yi 

/ 1^2 \ 0^2 V "2 

UJ 

2 

i ^ 

f Pl \ ai yrOLl 

G 2 J A_ 


L2 


conditioned on = ^, we have 
Pr (SIRo > P\uo G Uj,Yj = y) 


=Pr 

Mj-l 

= E 

n =0 




(^410 + ^410 +^4 


P 2 




n\ 


E 




ni,n2,n3j \Pj 


ni+n2 


^ E /,,10 


y,ic exp ( irkAC 


p. 

Pi 


X E/, 


plo exp -wkpo 


4|exp -/32/“^-^42 


M.-l 


= E 




n =0 


n\ 


E 


(^a)a=l 


n 

ni,n2,n3 


ni+n2 


where C^p\s^r) denotes the nth-order derivative of the 
Laplace transform Ci{s^r). 

Now, we calculate Ci{s^r) and C^p\s^r), respectively. 
First, ^^( 5 , can be calculated as follows: 


^ («’ r.'.io) ^ 




(’^s) / 


=E, 




exp(-s \D 


— CX\ 

1/01 gl,- 


{a) 

= Y^ic 


=E<^ 


ie^ic 

n E{gi,a exp(-s|Z)i,/or“'gi,£ 

1 


Ue^ic 

n 


= exp [ —27tP c 


1 + 5|Di^^o| 

(44,T2)Ai [ 

J r. 


1 - 


1 


1 + ^ 
I Coc I 


rdr j 

(31) 


where gp/ = 




, (a) is obtained by noting that 
(i G ^ic) are mutually independent, and (b) is obtained by 
using the probability generating functional of a PPP [10]. 
Further, by first letting s~^r = t and then = w, 

we obtain the result in ( p^ . 

Next, based on ( [sT] ) and utilizing Faa di Bruno’s formula 




































[11], 4ric(«> can be calculated as follows: ^ 0, i.e., 


nil 


nOO CO 

lim / y^Tj,Y,{n,y, U, Ti, T 2 , /y, {y)dy . (36) 


n (27rpe(C/,ri,T2)Ai [ 

a=l \ 


(’^a)a=l^-^rii 

rj,io 






•j.io V 1 “I" r“i 
ni!(-l) 


rdr 


We note that 

/ (1 — z)^ 

B (a, 6, 2 ;) =---h o ((1 — z)^) , as 2 : ^ 1 . (37) 


Then, we have 


(’^a)^=ieMni 


naii(a!)’"“’«a! 


ni / 

n 27rpe(C/,ri,T2)Air(a + l) / 

a=l \ 


pi—aai 


-dr 


\ o o 1 + cp 

\ o o 1 + cp 


ic0y 


1-^ 

a 

(c/5)“-^ 


o(p-») , (38) 

I - 1 ^ « / 2 ■ , (39) 

/i , s \a+i / V a ol + cpy a-- v / 

rj,ic (l+p2^j / ^ 

(32) where c G M+. Based on these two asymptotic expressions, 

where the integral can be solved using similar method obtaiid] 

as calculating £ 7 ^ ,^( 5 , r^pc, ^ipo). Similarly, we can cal- ^(ni) , x ^ 


ni! 


culate (s,r^po), ( 5 ,r^po), Ci.^{s,rj^ 2 ) and 

r^p). Finally, removing the conditions on ^ 

and after some algebraic manipulations, we can 
final result. 

D. Proof of Theorem 

Conditioned on = ^, we have 
1 - Pr (SIRo > /3|ixo e = y) 

/ 2 . 2 

=Pr 




KU^al 




+ o(^"^) , 


{s,y) = l3^ 


-(n2) 

n2 / 27r ' 


E 


n2! 


(Pa):iie^„ 


n:iiPa! 


d.oo^i.o 


< I3y°‘^ (rZlj,ic + -^Ij,io + 


Pi 


Pi 


X 


n 

a=l 


^Ai /pp^ p / 1 


a - ^ \ Pi 


= exp + ^Ij ,2 

~ (/jy«,p/p p p 

^ i + ^4io + -^4; 


4?:^ 4P = /3" 


n=M^ 
00 


E 


77 . 3 ! 


ns / 27r ' 

Oi2 ' 


= E 


i-l3y°‘^y 


n=Mj 
.(m) 


n! 


E 

(^a)a=i GAAn 


n 

ni,n2,n3 
(n2) 


Pi 


ni+n2 


n 

a=l 


(ga)a,= l^-^ri3 


ifA2 /P2A- 


n:ii9a! 


a - P VP 

^ a2 ^ 


<7a 


y' 


+ o(/3”^) 


(40) 


■o(/?"=) , 

(41) 


(42) 


4Pc (s>^jMC,ri,io) is,rj,io) 

( 33 ) 

Similar to the calculations in Appendix]^ after some algebraic 
manipulations and removing the condition on = ^, we have 


Moreover, utilizing dominated convergence theorem, we can 
show that 


poo ^ 

lim / y] Uy, (n, y, U, P) /y, {y)Ay 

n = M. 


pOO ^ 

= Y1 14 ,'^■-^3 P) 2/) P ^ 1 ) ^ 2 , P) fY, {y)dy . (43) 

0_ T\/r P 


n=M, 


l-5p[/,Ti,T2,P 


i 


3 

00 ^ 


where 


Hence, substituting ( [40| , •ED and ( [42| ) into and after 
some algebraic manipulations, we obtain Results 1), 2) and 
y] Tj,y, {n,y,U,TYT 2 ,P)fYyy)dy (34) 3) in Theorem^ To complete the proof, we now show that 

62 (f/, Ti, T 2 ) decreases with U. This can be proved by noting 
that i) 62 (P, Pi, P 2 ) is an increasing function of Pc (P, Pi, P 2 ), 
and ii) Pc ((7, Pi,P 2 ) decreases with U (which can be easily 
shown using ([29])). 


n=M- 


n 

ni,n2,n3 


Tj,Yj{n,y,U,Ti,T2,p) =— y] 

(^a ) a = l n 

^ y'^ L=/3y“3 ^4rio y) L=/3s/“). 4 

Now, we calculate the asymptotic outage probability when 


E. Proof of Lemma ^ 

1) Proof of Result 1): It can be easily shown that when 
(35) P ^ {0,1,..., A^i — A^ 2 }, we have Ni — U > N 2 . Thus, 

“/(a^) = o{g{x)) means lima^^o ^ = 0. 























d{U) = N 2 . When U G {A^i — N 2 , Ni — 1}, wc have 

Ni-U < N 2 . Thus, d{U) = Ni-U < N 2 . Hence, by @, 
we have = N 2 . 

2) Proof of Result 2): First, we note that Ti, T 2 ) 

can achieve the optimal order gain, i.e., Ti,T 2 ) G 

{0,1,...,A^i — A^ 2 }- Next, we compare the coefficients of 
which are achieved by different U G {0,l,...,A^i — 
A^ 2 }- We consider two cases, i) When U < Ni — 
N 2 , as b 2 {U,Ti,T 2 ) decreases with U, the coefficients 
are 62 (A^i - A ^2 - 1 , Ti, T 2 ) < 62 (A^i - Ar 2 - 2, Ti, T 2 ) < 

... < 62 ( 0 ,Ti,T 2 ). ii) When U = A^i — A^ 2 , the coefficient 
of (3^- is 61 (ATi - N 2 , Ti, Ts) + b 2 (fVi - N 2 , Ti, T 2 ). 

Based on the discussions above and ( |^ , we 
know that the exact value of Ti,T 2 ) depends 

on whether 62 (A^i — A ^2 — 1 , Ti, T 2 ) is smaller than 
61 (A/^i — A ^2 5 Ti,T 2 ) + ^2 (-^1 ~-^ 2 : Ti, T 2 ) or not. 
In particular, when 62 (^1 — ^2 — 1, Ti, T 2 ) < 

61 (A^i — A^ 2 ,Ti,T 2 ) + 62 (A^i — A/' 2 ,Ti,T 2 ), we 

have Ti, T 2 ) = Ni — N 2 — 1 ; otherwise, 

f/*(ATi,T2)=7Vi-7V2. 


